Observation of the distribution of molecular spin states 
by resonant quantum tunneling of the magnetization 
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Below 360 mK, Feg magnetic molecular clusters are in the pure quantum relaxation regime and 
we show that the predicted "square-root time" relaxation is obeyed, allowing us to develop a new 
method for watching the evolution of the distribution of molecular spin states in the sample. We 
measure as a function of applied field H the statistical distribution P{(,h) of magnetic energy bias 
acting on the molecules. Tunneling initially causes rapid transitions of molecules, thereby "digging 
a hole" in P{^h) (around the resonant condition ~ 0). For small initial magnetization values, 
the hole width shows an intrinsic broadening which may be due to nuclear spins. 



PACS numbers: 75.45.+j, 75.60Ej 

Strong evidence now exists for thermally-activated 
quantum tunneling of the magnetization (QTM) in mag- 
netic molecules such as Mni2ac and Fes Crystals of 
these materials can be thought of as ensembles of identi- 
cal, iso-oriented nanomagnets of net spin 5 = 10 for both 
Mni2ac and Feg, and with a strong Ising-like anisotropy. 
The energy barrier between the two lowest lying spin 
states with Sz — ±10 is about 60 K for Mni2ac and 
25 K for Feg Theoretical discussion of thermally- 
activated QTM assumes that thermal processes (princi- 
pally phonons) promote the molecules up to high lev- 
els, not far below the top of the energy barrier, and the 
molecules then tunnel inelastically to the other side. The 
transitions are therefore almost entirely accomplished via 
thermal excitations. 

At temperatures below 360 mK, Feg molecular clus- 
ters display a clear crossover from thermally activated re- 
laxation to a temperature independent quantum regime, 
with a pronounced resonance structure of the relaxation 
time as a function of the external field This can 
be seen for example by hysteresis loop measurements 
(Fig. |l|) . In this regime only the two lowest levels of each 
molecule are occupied, and only "pure" quantum tunnel- 
ing through the anisotropy barrier can cause direct tran- 
sitions between these two states. It was surprising how- 
ever that the observed relaxation of the magnetization in 
the quantum regime was found to be non-exponential and 
the resonance width orders of magnitude too large . 
The key to understanding this seemingly anomalous be- 
havior now appears to involve the ubiquitous hypcrfinc 
fields as well as the (inevitable) evolving distribution of 
the weak dipole fields of the nanomagnets themselves . 

In this letter, we focus on the low temperature and low 
field limits, where phonon-mediated relaxation is astro- 
nomically long and can be neglected. In this limit, the 



Sz = ±10 spin states are coupled by a tunneling matrix 
element Atunnci which is estimated to be about 10~* K 
1^ . In order to tunnel between these states, the magnetic 
energy bias = QIi-bSH due to the local magnetic field 
H on a molecule must be smaller than Atunnei implying 
a local field smaller than 10^^ T for Feg clusters. Since 
the typical intermolecular dipole fields are of the order of 
0.05 T, it seems at first that almost all molecules should 
be blocked from tunneling by a very large energy bias. 
Prokofev and Stamp have proposed a solution to this 
dilemna by assuming that fast dynamic nuclear fluctua- 
tions broaden the resonance, and the gradual adjustment 
of the dipole fields in the sample caused by the tunneling, 
brings other molecules into resonance and allows contin- 
uous relaxation A crucial prediction of the theory is 
that at a given longitudinal applied field H, the magne- 
tization should relax at short times with a square-root 
time dependence: 



M{H, t) = Mi„ + (Meq(i7) - Min)y^r,qrt(CH)i (1) 

Here Afjn is the initial magnetization at time t — Q (i.e. 
after a rapid field change), and Mcq{H) is the equilib- 
rium magnetization. The rate function rsqrt(^ff) is pro- 
portional to the normalized distribution P{^h) of energy 
bias in the sample: 



rsqit(fff) 



(2) 



where h is Planck's constant and c is a constant of the 
order of unity which depends on the sample shape. If 
these simple relations are true, then measurements of 
the short time relaxation as a function of the applied 
field H gives experimentalist a powerful new method to 
directly observe the distribution P{£^h)- Indeed the pre- 
dicted y/t relaxation (Eq. (|l|)) has been seen in prelimi- 
nary experiments on fully saturated Feg crystals . We 
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show here that it is accurately obeyed for saturated and 
non-saturated samples (Fig. |^) and we find that a re- 
markable structure emerges in P{^h) as presented in the 
following. 

In order to carefully study P{£,h) and its evolution as 
the sample relaxes, we have developed a unique magne- 
tometer consisting of an array of micro-SQUIDs @jll) 
on which we placed a single crystal of Feg molecular clus- 
ters. The SQUIDs measure the magnetic field induced by 
the magnetization of the crystal (see inset of Fig. |l|) . The 
advantage of this magnetometer lies mainly in its high 
sensitivity and fast response, allowing short-time mea- 
surements down to 1 ms. Furthermore the magnetic field 
can be changed rapidly and along any direction. 

Figure 1^ shows a typical set of relaxation curves plotted 
against the square-root of time. However instead of sat- 
urating the sample before each relaxation measurement 
so that the initial magnetization Afjn = Mg as described 
in 1^ , these measurements were made by rapidly quench- 
ing the sample from 2 K in zero field (ZFC), i.e. for an 
initial magnetization Mjn = 0. The quench takes approx- 
imately one second and thus the sample does not have 
time to relax, either by thermal activation or by quantum 
transitions, so that the high temperature "thermal equi- 
librium" spin distribution is effectively frozen in. Once 
the temperature is stable (in this case 40 mK) a mea- 
suring field is applied, the timer is set to i = 0, and the 
relaxation of the magnetization is recorded as a func- 
tion of time. The entire procedure was repeated for each 
measuring field shown in Fig. |^. As can be seen for short 
times t < 100 s the square root relaxation is well obeyed. 
Note that all curves extrapolate back to Af = at i = 0. 
A fit of the data to Eq. (|l|) determines Fsqrt- 

A plot of Fsqrt vs. H is shown in Fig. || for the zero 
field cooled data, as well as distributions for three other 
values of the initial magnetization which were obtained 
by quenching in small fixed fields (field cooled FC mag- 
netization). The distribution or an initial magnetization 
close to the saturation value is clearly the most narrow re- 
flecting the high degree of order starting from this state. 
The distributions become more broad as the initial mag- 
netization becomes smaller reflecting the random fraction 
of reversed spins. The small satellite bumps are due to 
flipped nearest neighbor spins on the tri-clinic lattice as 
seen in computer simulations [^,^ . 

We can exploit this technique of measuring P{^h) in 
order to observe the evolution of molecular states in the 
sample during relaxation by quantum tunneling of the 
magnetization. We flrst fleld cooled the sample (ther- 
mally anneal) as described above in order to obtain the 
desired initial magnetization state. Then after applying 
a field Ht, we let the sample relax for a time tt, which 
we call "tunneling field" and "tunneling time" respec- 
tively. During the tunneling time, a small fraction of the 
molecular spins tunnel and reverse their direction. Fi- 
nally, we applied a small measuring field and record the 



short time relaxation which again can be fit to a square 
root law (Eq. (|^)) yielding Fsqrt- The entire procedure is 
repeated many times for other measuring fields in order 
to probe the distribution as a function of field i?, and 
thus we obtain the distribution P{£^h i Ht,it) which we 
call a "tunneling distribution" . 

Figure ^ shows tunneling distributions for field Ht — 
and for tunneling times between 1 and 250 s for the 
case that the initial magnetization starts from the fully 
saturated state. Note the rapid depletion of molecular 
spin states around the resonant field Ht = and how 
quickly the depletion depth and width increase with tun- 
neling time. In effect, a "hole is dug" into the distribu- 
tion function around Ht- The hole arises because only 
spins in resonance can tunnel. The hole is spread out 
because as the sample relaxes, the internal fields in the 
sample change such that spins which were close to the 
resonance condition may actually be brought into reso- 
nance. Notice however that "wings" are created on each 
side of the hole because other spins are pushed further 
away from resonance. These features are in good agree- 
ment with Monte Carlo simulations of the relaxation for 
non-spherical samples [p|,p|, p"2|jT^ ] . 

In Fig. |^(b) we see the extraordinary effect of sample 
annealing (i.e. for small values of the initial magneti- 
zation Mjn) on the evolution of P{^h , Ht,tt) with time. 
Now the depletion proceeds over an extremely narrow 
bias range. This is virtually incontestable experimental 
proof that we are seeing tunneling relaxation. The nar- 
rowing of the hole is because in the annealed sample fur- 
ther incremental relaxation hardly changes the internal 
demagnetization field. Notice that the initial line shape 
of P{^h) is very accurately fit to a Gaussian for the an- 
nealed samples, exactly as predicted for the dipole field 
distribution of a dense set of randomly oriented spins 

0- . . . 

Further investigation of the effect of sample annealing 

led us to another remarkable discovery (Fig. ^). Pro- 
gressive annealing such that |Min| < |0.5Ms|, eventually 
leads to a hole linewidth which at short times is inde- 
pendent of further annealing, and has a half linewidth of 
0.8 mT. It is interesting that such an intrinsic linewidth 
was predicted by Prokof'ev-Stamp ||]. It is claimed to 
come from the nuclear spins which would give rise to a 
linewidth of roughly the same order (although only 2% 
of natural iron has a nuclear moment, there are other nu- 
clei in the clusters that can contribute to the hyperfine 
flelds, i.e. more than 100 hydrogen, 18 nitrogen and 8 
bromine atoms!). We notice that any intrinsic linewidth 
due to the tunneling matrix element itself is 5 orders 
of magnitude smaller, and would be quite unobservable. 
According to Eq. (H), the ratio Fsqrt /-Ed (where Ejj is 
the Gaussian half-width of P{£,h) for strongly annealed 
samples) should be a constant, and thus allows us to esti- 
mate Atunnci from our relaxation measurements. We flnd 
that it is indeed a constant (even though Ejj and Fgqrt 
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vary with Min), and we extract Atunnei « 5 x 10~* K for 
|Min| < |0.5Ms|, assuming c = 1. This agrees well with 
the expected value |^,^. 

In conclusion, we have developed a new measurement 
technique yielding P{^h) which is related to the inter- 
nal dipole field distributions always present in crystals of 
molecular clusters. The distribution evolves during relax- 
ation by tunneling in a non-trivial way, and can be mon- 
itored by our technique, revealing the details of how the 
tunneling is proceeding in the sample, which molecules 
are tunneling, and how the time- varying internal fields in- 
fluence the relaxation. The shape of the hole for thermal 
annealed distributions indicates a fast dynamic relax- 
ation over a field range of 0.8 niT which could correspond 
in the Prokof'ev-Stamp theory to the nuclear linewidth 
^0 • Although this is only indirect evidence of the nuclear 
mechanism, it is hard for us to see what else could be 
operating at these temperatures. Our evidence for the 
role of the dipole interactions is on the other hand very 
direct, and in good agreement with Monte Carlo simu- 
lations. We believe that our technique should work for 
other multi-particle spin systems in the quantum regime 
(like quantum spin glasses ) , and could give quite new 
information on the non-ergodic relaxation behavior typ- 
ical of these systems. 
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FIG. 1. Magnetic hysteresis curves for a crystal of Feg 
molecular clusters in the quantum regime (40 mK) for three 
field sweeping rates. Resonant tunneling is evidenced by six 
equally separated steps. Note the strong dependence on the 
field sweeping rate. The measurements were made using an 
array of micro-SQUIDs (see inset) [ p"o|jll| ]. The high sensitiv- 
ity of this magnetometer allows us to study single crystals of 
the order of 10 to 500 fim which are placed directly on the 
array. 

FIG. 2. Typical square root of time relaxation curves for 
an Feg crystal measured at 40 mK. For each curve, the sam- 
ple was first thermally annealed at if = (ZFC), then a 
small field was applied and the relaxation of magnetization 
was measured during 1000 s. The slope of the lines gives 
Fsqrt when plotted against the square-root of t as shown. 
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FIG. 3. Field dependence of the short time square root re- 
laxation rates raqrt(5jf) for three different values of the initial 
magnetization Min. According to Eq. the curves are pro- 
portional to the distribution P{^h) of magnetic energy bias 
£,H due to local dipole field distributions in the sample. Note 
the logarithmic scale for Fsqrt. The peaked distribution la- 
beled Min = —0.998A4s was obtained by saturating the sam- 
ple, whereas the other distributions were obtained by thermal 
annealing (FC or ZFC). Mn = 0.870Ms is distorted by near- 
est neighbor lattice effects. 
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FIG. 4. Tunneling distributions: the field dependence of 
the short time square root relaxation rates rsqrt(^H) are pre- 
sented on a logarithmic scale showing the depletion of the 
molecular spin states by quantum tunneling at _ff = for var- 
ious tunneling times. 

(a) Tunneling distributions for the initial magnetization start- 
ing from saturation Afin = Ms. For each point, the sample 
was first saturated in a field of -1.4 T and at a temperature of 
about 2 K and then rapidly cooled to 40 mK. After applying 
a "tunneling field" of Ht — 0, we let the sample relax for a 
"tunneling time" tt (after tt = 2, 10, 50 and 250 s, the re- 
versed fraction of magnetization is 0.008, 0.030, 0.077, 0.159 
of Ms , respectively) . Finally, we applied a small field to mea- 
sure the short time relaxation which could be fit accurately 
to a square root law yielding Fsqrt. Because Fsqrt is propor- 
tional to the spins which are still free to tunnel, one obtains 
the distribution P{£,h , Ht,tt). (b) Tunneling distributions as 
in (a), but now for each point, the sample was first annealed 
(FC) to a value Mn = -0.2Ms. After U = 5, 10, 20 and 40 s, 
the reversed fraction of the magnetization is 0.0012, 0.0020, 
0.0032, 0.0049 of Ms, respectively. At the resonance, the de- 
pletion develops very rapidly with elapsed time, even though 
the total magnetization and the internal demagnetization field 
hardly change during this time. Notice that parabolic shape 
of Fsqrt (^h) shows it is accurately Gaussian, with a half- width 
of 0.030 T. This Gaussian profile is found for |A/in| < \0.5Ms\ 
(but with a half- width Ed depending on A/in). 

FIG. 5. Detail showing the depletion of molecular spin 
sates by quantum tunneling for different values of annealing 
(FC). The difference between the initial thermally quenched 
distribution (i.e. tt = 0) and the tunneling distribution ob- 
tained after allowing the sample to relax for = 16 s are 
plotted. For initial FC magnetization values close to satu- 
ration, the depletion is larger and asymmetric whereas it is 
narrow and symmetric for initial magnetization values which 
are smaller than |0.5Afs|. Note that in the latter case the hole 
is independent of initial magnetization, yielding an intrinsic 
hole width of about 0.8 mT (see also fig. Bj) . 
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